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1. Introduction 

Numerical simulations of fermion systems have to deal with the Pauli principle which enforces 
a completely anti-symmetric wave function for fermions. Obviously, this requirement makes 
fermion systems extremely non-local. Changing the system at a single point affects all degrees 
of freedom. In more technical terms, an ad-hoc local change may lead to a completely different 
value of the fermion determinant. Thus intricate methods such as the hybrid Monte Carlo algorithm 
were developed. Nevertheless, also these methods have to deal with the non-locality which makes 
the simulation of fermions several orders of magnitude more expensive than bosonic systems. 

Circumventing the non-locality problem altogether is certainly an extremely appealing idea. 
A prominent example of a breakthrough in this direction is the Meron Cluster Algorithm proposed 
in [|l|], which allows for highly effective Monte Carlo simulation for certain classes of fermionic 
models. 

An alternative approach is a formal solution of the Grassmann path integral for fermions which 
represents the partition function of the system as a model of closed fermion loops. While for a long 
time this technique has been known to work well in the strong coupling limit, only recently [^, ||] 
loop representations were found for two-dimensional lattice field theories at arbitrary coupling, in 
particular the Gross-Neveu model [Q]. 

In a recent publication |Q] it was shown that the loop representation allows for an efficient 
and considerably cheaper simulation than traditional methods. In a subsequent paper [^] Wolff has 
rederived the loop representation by decomposing 2-d Dirac fermions into Majorana components 
and demonstrated that the loop formulation can be recast as a spin system where a cluster algorithm 
boosts the efficiency of a numerical simulation further. 

In this contribution we review the loop representation of the lattice Gross-Neveu model and its 
use for a numerical simulation. We furthermore discuss the loop representation of the Schwinger 
Model [^, i.e., QED in two dimensions, an example which illustrates the limitations of the loop 
approach for a use in numerical simulations. 

2. Loop representation of the lattice Gross-Neveu model 

We consider the lattice Gross-Neveu model with A'^ flavors of Wilson fermions. The lattice action 
for the fermions is given by (we set the lattice spacing to a = 1) 

D{x,y) = [2 + m + (p(x)] 5,-,^ - £ li^^^^^^ . (2.1) 

Ai=±l ^ 

The sum runs over the 2-d lattice A. In two dimensions the /-matrices may be chosen as the Pauli 
matrices, = o^. The spinors 1/7 and i/A are vectors of N 2-d spinors, and we use vector/matrix 
notation for both the spinor and flavor indices. Through the Dirac operator D{x,y) all flavors couple 
in the same way to the real scalar field (p, which has the action 

SsM = ^"L ^^^f ■ (2.2) 
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When integrating out the scalar field the 4-fermi interaction 



A-eA 



(2.3) 



is induced. The partition function of the model is given by 

Z= /p]^/<p(;c)^/1^7(;c)^/v^(;c)^'-^sW-^^[^■V'''Pl = jY\d(p{x)e-^^^'^Ua{D[(p]Y , (2.4) 

X X 

where in the second step the fermions were integrated out giving rise to a remaining path integral 
over the scalar field with the fermion determinant raised to the power N as integrand. 

The last expression is a possible starting point for identifying the loop representation. The 
Dirac operator D may be rewritten as 



D{x,y) = [2 + m + (p{x)] 5,^y - H{x,y) 
1 



(2.5) 



H{x,y) 



y ^ 



2+m+ (p{x) 

where we have combined all nearest neighbor terms in the hopping matrix H. Inserting the repre- 



sentation ( |2.5[ ) into the partition function ( |2.4| ) one finds 

Z = /"n^/(^»(x)e-^^['''ln(2 + "J + <PW)^'^det(l-//[(p])' 



Yld(l>{x)e-^'^f^ Yl{2 + m + (pix)f^ ( exp ( - £ -Tr 



n=l 



H" 



(2.6) 



where we have used the formula det[l —H]= exp(Tr ln[l —H]) for the determinant and expanded 
the logarithm. 

The expression ( ^^ is the well known hopping expansion. At this point the loops are al- 
ready evident: The hopping matrix // is a matrix which describes hopping between neighboring 



lattice points. Consequently the power H" in ( |2.6[ ) corresponds to a chain of n subsequent steps. 
When taking the trace only closed chains, i.e., loops survive. Such an expansion holds in arbitrary 
dimensions and for different types of bosonic fields, scalar, as well as gauge fields. 

The crucial step, however, is that the traces Tr [H"] in ( ^^ can be evaluated only in special 
cases. In addition to the space-time indices, this trace is over Dirac and for non-abelian gauge 
theories also over the color indices. For the latter a simple closed form is probably not realistic. 
Concerning the Dirac indices, in two dimensions it is possible [§] to find a closed form for the trace 



over the matrices [1 it 7^]/2 which enter the hopping matrix (p3|). Thus for the case of non-abelian 



interactions in two dimensions the exponent in ( [2.6[ ) can be computed in closed form. 

As discussed, for special cases (scalar or abelian bosonic fields in 2-d) the coefficients for the 



individual loops in the exponent of ( |2.6| ) can be computed analytically. The final step is to bring 
the loops down from the exponent. Here two different approaches were followed in ^ and [^. In 
the former case the final expression for the loop representation was obtained by comparing the 2-d 
Wilson fermions to the hopping expansion of a 8-vertex model. In the latter case a direct evalua- 
tion of the exponential of the sum over loops was performed. Again we remark that there is also 
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the direct identification of the loop representation through the expUcit solution of the Grassmann 
integral for the Majorana components 

Once the determinant is given as a sum over loops with known coefficients (not as the expo- 
nential of a sum over loops) there is only the path integration over the bosonic variables attached to 
the loops left to be done. For the case of the scalar fields which give rise to the 4-fermi interaction, 
the path integration is trivial, since at each lattice point only moments of the Gaussian distribution 
need to be computed. In this way the partition function of the A^-fiavor lattice Gross-Neveu model 
is found to be a model of 2N self-avoiding loops. For the case of general N we refer the reader to 
[^, and here quote the result forN =\, which is the case that was used in the numerical simulations 
[Q, ^]. The partition function reads 



r.b 



The sum runs over two sets of loops which we refer to as red (r) and blue {b). For a given color 
the loops are self avoiding, i.e., they cannot cross or touch each other, while loops of different 
may do so. In Eq. ( [2.7[ ) c{r,b) is the total number of corners for both, red and blue loops. Thus 
every comer contributes a factor of l/\/2 to the weight of a configuration. Furthermore, n\{r,b) 
is the number of lattice sites which are singly occupied by either r or b and n2{r,b) is the number 
of doubly occupied sites, i.e., sites which are visited by both, a red and a blue loop. The weight 
factors /i and /2 are simple functions, related to the mass m and the coupling g through 

_ 2+m _ 1 

^1 = (2 + m)2 + g ' = (2 + m)2 + g " ^^"^^ 



The mapping (2.7), (2.8) is exact in the thermodynamic limit. For finite volume different types of 
boundary conditions in the two representations lead to finite size effects: In the loop representa- 
tion we need to have closed loops and in a finite volume the loops can wind around the periodic 
boundary. The loop configurations fall into three equivalence classes, C^"^ jC^" depending on 
the numbers of red and blue non-trivially winding loops (see also |^): C'^'^ (even-even): The to- 
tal number of windings for both, red and blue loops is even for both directions. C^" (even-odd): 
One of the colors has an odd number of windings for one of the directions. C"" (odd-odd): Both 
colors have an odd number of windings in one of the directions. These equivalence classes can- 
not be linked in a simple way to the boundary conditions in the standard representation which we 
discussed above. However, in it was shown that the boundary effects vanish as 1/ \/V, with V 
denoting the volume. The representation in terms of the Ising spin variables ||^ solves the boundary 
condition problem completely, and the partition functions of the original fermionic- and the spin 
representation are identical also on finite volumes. 



3. Numerical simulation 



For the numerical simulation of the loop representation of the Gross-Neveu model we use a 
local Metropolis update. Red and blue loops are updated alternately by performing a full sweep 
through the lattice for only one color and meanwhile treating the other as a constant background 
field. During one sweep all plaquettes are visited once. A trial configuration is offered by inverting 
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the 4 links of the current plaquette of the active color With such an offer, we guarantee that the 
loops stay closed, or new loops are created if all links of the plaquette were empty before. In 
case that the self-avoidance condition is violated, the proposal is rejected. Otherwise the new 
configuration is accepted with the Metropolis probability 

/f"'/f"'. (3.1) 

Ac is the difference of the number of comers, Awi and A«2 are the differences in the occupation 
numbers. Fig. |l| shows snapshots of typical loop configurations in the numerical simulation. 

Particularly simple observables are derivatives of the free energy F = — \nZ. These expres- 
sions can be written as moments of occupation numbers. To be more explicit we discuss the chiral 
condensate % and its susceptibility C^. The conventional definitions are 

x = yl^mxmx)) = -y^. (3.2) 



dm 



Cx = ^- (3-3) 



In terms of loop variables these expressions read 

^ = -V7 [/2(«i)+2/f(«o)] , (3.4) 



+ 2/iV2((«o + «i-(«o + «i))'> - (4/f-2/iV2)(«o)-/2'(«i)] , (3.5) 
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Figure 2: L.h.s.: The chiral condensate x for g = as a function of m for 2 different lattice sizes. We 
compare the simulation in the loop representation (symbols with error bars) to the exact result from Fourier 
transformation (curves). R.h.s.: Same as on the l.h.s., now for the chiral susceptibility €■),. 



where the number hq is the total number of empty lattice sites. These representations were obtained 
by differentiating the partition function (2.7). 

In Fig. |2| we compare the loop results at ^ = (symbols) with those from Fourier transfor- 
mation (curves) which is applicable in standard representation for the special case of vanishing 
coupling. Two volumes are used, the relatively small lattice 32 x 32 and a substantially larger lat- 
tice of 512 X 512. For each value of m typically 10000 sweeps were used to equilibrate the system 
and then about 50000 measurements were performed for our observables. These measurements 
have been separated by 10 sweeps for each color to guarantee sufficient decorrelation. For the cal- 
culation of the observables we used Eqs. (3^), (3^), and the statistical error was computed with 
the jackknife method. For the larger lattice we find almost perfect agreement of the results from 
the loop representation with the analytic results. In [^] we have presented the results from the loop 
simulation for several values of ^ 7^ and we compared these to the outcome of a simulation with 
traditional techniques. Concerning the performance we find that with comparable amounts of CPU 
time with the loop representation we can work on volumes which are two orders of magnitude 
larger than those available with traditional techniques. The cluster approach of iQ] enhances this 
performance further. 



4. Perspectives and limitations of the loop representation 

Having addressed the merits of the loop representation for the Gross-Neveu model, we would 
like to comment on possible extensions of the loop method, but also discuss the points where we 
see limits of the method. 

We begin this discussion with stressing that, although we so far restricted our numerical sim- 
ulations to only two flavors, a generalization to the Gross-Neveu model with an arbitrary number 
of flavors is straightforward with the loop formula given in [Q]. 

Concerning models with relativistic fermions which are coupled via 4-fermi interactions in 
higher dimensions, one could try to repeat the strategy that led to the loop representation for the 
2-d models. An essential step in the identification of the loop formalism was the closed result for 
the traces of the 7-matrices. While this is a relatively simple problem in 2-d, the corresponding 
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Structures in 4-d are considerably more involved. For an attempt to find such a closed formula in 
four dimensions see, e.g., [Il0|]. 

Interesting might also be the case of non-relativistic fermions in 2+1 dimensions with 4-fermi 



interaction. For some of these systems a relation to 3-d spin models is known [[111, [12|] which might 
be useful for a numerical simulation. 

We finally comment on the applicability of the loop approach to lattice gauge theories beyond 
the strong coupling limit. We have already mentioned, that we judge the case of non-abelian gauge 
fields as an elusive goal, due to the non-commutativity of the link variables. For abelian gauge 



fields the situation is simpler and in 2-dimensions the loops in (2.6) can again be computed in 
closed form. The resulting loop representation for the lattice Schwinger Model ||7|] is of a different 
type, however. Since gauge fields are oriented quantities, one has to work with oriented loops. 



while the loops for a scalar interaction are non-oriented (see Eq. (|2.7[)). One finds that reverting 
the orientation of a loop corresponds to complex conjugation of its contribution. This implies, that 
certain cancellations among loops, which simplify the scalar case, are no longer possible 0. The 
loops for the Schwinger Model turn out to be self-intersecting and an extra minus sign appears 



for each intersection. In a numerical simulation [ ]13[ ] it was found that the resulting fermion sign 



problem limits the size of the accessible volumes. At the moment it is unclear whether this is a 
fundamental obstacle or if this problem can be overcome by different techniques. 
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